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Abstract 

We introduce a simple yet powerful calculational tool useful in calculating averages of ratios and products 
of characteristic polynomials. The method is based on Dyson Brownian motion and Grassmann integration 
formula for determinants. It is intended as an alternative to other RMT techniques applicable to general 
gaussian measures. Resulting formulas are exact for finite matrix size N and form integral representations 
convenient for large N asymptotics. Quantities obtained by the method can be interpreted as averages over 
matrix models with an external source. We provide several explicit and novel calculations showing a range 
of applications. 


1. Introduction 


One of the strengths of Random Matrix Theory lies in the abundance of powerful calculational tools, with 
the method of orthogonal polynomials [i| , supersymmetric techniques Q and free probability Q among many 
others. This paper attempts to enlarge this family with a convenient technique we call the diffusion method. 
It serves as a fast and powerful framework useful when dealing with powers and ratios of characteristic 
polynomials averaged over external source gaussian measures. It began as a byproduct of QGD considerations 
made several years ago and was thereafter successfully applied to hermitian, Wishart and chiral models 
0 S 0 • The method uses a Dyson-like picture of dynamical matrices to achieve this goal. 

Studying characteristic polynomials in the RMT community is now a prolific topic with many branches, 
its root however can be traced back to a remarkable relation connecting them to orthogonal polynomials 


(det(2 - X))^ = 


Such objects were considered in [8| and in many areas of application such as zeroes of Riemann C, function [9j, 
eigenvalue statistics in quantum chaotic systems and matrix models of QGD [h|. Moreov er, p roducts 
and ratios of characteristic polynomials reveal rich mathematical structures in both hermitian 1^ [13| and 
non-hermitian IJ, Il5l . Il6l | models. 


At the core of the method lies a seminal work of Dyson [17| who observed that a static matrix model 
can be interpreted as a dynamical system. He showed that the joint probability density function for N 
eigenvalues behaves exactly like a statistical system of N ’’particles” interacting via logarithmic potential. 
Such particle system undergoes a Dysonian Brownian motion defined by the Langevin equation of the form; 


d\i — 


1 




-dt. 


This interacting Brownian motion is induced by a gaussian diffusion applied independently to matrix entries. 
Even though it lacks the harmonic potential usually present to attain the stationary the r —>■ oo limit, we 
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show how to include this canonical approach of Dyson in the present framework. Studying RMT from the 
Brownian motion’s point of view attracted attention of physicists 18, El and mathematicians (HElj alike. 

Plan of the paper is as follows. In Sec. 2 we discuss the framework of the method in full generality and 
comment on its characteristics. In Sec. 3 we comment on the relation to standard random matrix models 
and on how to include the harmonic potential. In Sec. 4 we calculate five examples with special attention 
given to non-hermitian models. We show how to arrive at the formula for the averaged ratio of characteristic 
polynomials wrt. (3 = 2 complex hermitian models and derive a novel duality-type equation for averaged 
products of characteristic polynomials in the complex (3 = 2 non-hermitian ensemble. Furthermore, we 
compute a new integral representation of the averaged characteristic polynomial in a correlated variation of 
complex (3 = 2 non-hermitian ensemble, compute the same object for a multiplication of two non-hermitian 
matrices and study real/complex non-hermitian ensemble crossover. The examples provided in the last 
section are mostly novel results. 


2. Diffusion method 

We introduce entrywise diffusive dynamics to M - a, N x N matrix of interest with prescribed symmetries. 
Well-suited formalism for our purpose is the multidimensional heat equation: 

drPiM,T) = ^AmP{M,t), (1) 

where P{M,t) is the joint probability density function. Am denotes a Laplace operator over indepedent 
degrees of freedom of M and the constant 1/iV is a convention. For example, if M = Afl, the hermitian 
laplacian is flat and reads Am = 5 ^xu + \ ^ class of objects to which the method 

applies are the ratios and products of characteristic polynomials denoted in general by D{Z,M). We are 
ultimately interested in formulas for the average 

'd;{Z) := {D{Z,Mr))M^ = j dMP{M,T)D{Z,PI), (2) 

where {)j^ is a normalized averaging over the dynamical matrices Mr- Second equality is a consequence of 
the dehntion of jPDF P(X,t) = {S{Mr — 

To proceed, we extend D{Z,M) 'D(Z, M; A) by introducing parameter-like variables A such that 
limA-j-Ao D(Z, M; A) = D{Z,M) with Aq = 0 in most cases. This procedure in practice means we rewrite 
D as a block determinant and introduce off-diagonal terms accordingly. The sole purpose of this auxiliary 
construction is the search of diffusive equations in A-parameter space for the averaged deformed quantity 
D;(Z;A) = (D(Z,A/,;A))^^. 

Constructing useful deformation parameters A is not difficult however it varies between problems and is 
found in a case-by-case manner. Usefulness is determined solely by its ability to produce a closed diffusive 
type equation for D,-. In simplest cases introduction of additional parameters is superfluous since already 
the arguments Z of charactetistic polynomials suffice in deriving the diffusion equation. 

In order to find the final equation we consider a time derivative of T>r ' 

dr^ =jjJ dMAMP{M, r)D(Z, M;A) = ^J dMP{M, t)Am'D{Z, M; A), (3) 

where we used the equation ([I]) and integrated by parts to move the differential operator to D. Remaining 
task is to find Aa such that the condition 

AmD(Z, M- A) = AaD(Z, M; A) (4) 

is satished. We then write the diffusive equation as 

dr^iZ- A) = 1 AaD;(Z; A). (5) 
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As can be seen from the condition Q , gaussian laplacians on M manifold transform into Gaussian laplacian 
operators on A space but, at the same time, we observe a drastic decrease in the number of variables. 

Practical calculations solving use the Grassmann/complex representation of determinants or fermion- 
ization technique 


detM-y'de-^''^^’', detM-i - (6) 

where the constants are not important for our purpose. In this representation the derivatives are exception¬ 
ally simple to compute. This is also the reason why it is desirable of A deformations to retain determinantal 
structure. 

As a last step, we solve the equation ([S]) and set auxiliary parameters back to their initial values A —> Aq 
to recover the object of interest Dr- 


2.1. Comments and resume 

The method is applicable to any type of entrywise diffusion © which is perfect for models with nontrivial 
variance structure. An instance of such model was considered in Example 3. The diffusion equation in 
parameter space © has in general low dimensionality when compared to the matrix size, this duality 
phenomenon is widely observed in other methods. The resulting formulas, being the solution to a differential 
equation, are always dependent on the initial condition - the starting matrix Mg. For standard matrix models 
this is translated into a jPDF with an external source 22|, an observation elucidated in Sec. 3. Final formulas 
can be viewed also as integral representations convenient for large N analysis. 

A general way to proceed follows these subsequent steps: 

1 . introduce entrywise diffusion of choice o, 

2. define object D (i.e. product and ratios of determinants) and form a A parameter extension D, 

3. infer a diffusion equation in the A space for the averaged quantity D with the help of Grassmann/complex 
representation of D and condition (|T|) , 

4. solve the equation ® and set A parameters to its undeformed values Ag to recover the object of 
interest. 


3. Relation to standard random matrix models 

In previous section we have discussed a general framework in the diffusive language, here we comment 
on how to connect this approach to standard random matrix measures. Entrywise diffusion © is, as a 
differential equation, implicitly reinforced with initial conditions which we consider to be of delta function 
type P{M, r —>■ 0) = S{M — Mg). The matrix at time r is equal to 

Mr = Mg -h yPPM, (7) 

where At is a matrix chosen randomly from jPDF PMo=o(At, r = 1). By inverting this equation we conclude 
that 


P{M)dM ^ P dM{T), (8) 

where on the RHS the time r serves as a parameter and can be set to r = 1. A family of matrix models of 
the form P(M — Mg) are called external source or shifted mean models [2^. In present context, averaging 
over dynamical matrices with time r is equivalent to matrix models of variance proportional to r and with 
an external source: 


,T = 1 


X{M). 


dMP(M,T)X(M) 


dMP 


M - Mo 


( 9 ) 





The proportionality constant depends on the form of P. In particular, if M = = H is complex hermitian, 

the solution to (P) reads 


P{H,T) = Cexp(^-^TY{H-Hofy ( 10 ) 

with C = ■ In the case oi M = X complex without any additional symmetries we calculate that 

P{X,t) = C'exp (^-^Tr(X - Xo)HX - , (11) 

with C = {^) so that J PdX = 1. We observe how the time r serves now as scale of the variance and 
the external source being the initial matrix Xq{Hq). 

We now address the lack of harmonic potential in the entrywise diffusion (ED present in most treatments 
of Dyson Brownian motion. It serves as an effective confining well giving a non-zero stationary r —>• oo limit. 
Its presence can be however mimicked by a suitable reparametrization (^ . Namely, by setting the source 
Mq —>■ Moe“'^ and time variable r —|— we effectively retrieve a diffusion in the harmonic potential. 
As r —>■ oo, the initial matrix information is lost as the diffusion fills the whole confining potential which 
was already observed in Dyson’s paper [T^ . 


4. Examples 


This section is devoted to several interesting examples. It serves also as a tour-de-force showing the 
framework at work to calculate some new results and compare to known ones. The majority of them deal 
with the non-hermitian ensembles. 

Example 1 is devoted to the probably most thoroughly studied random matrix ensemble of complex 
P = 2 hermitian matrices. We show the applicability of our method to the averaged ratio of determinants, 
obtain an integral representation for any external source Hq and show how it reduces to known results [l^ 
for Ho —>• 0. 

Example 2 elucidates on a certain duality-type formula for complex /3 = 2 non-hermitian ensembles, a 
result which continues the successful programme of dualities obtained in hermitian 2J,[2^ and non-hermitian 
0 ensembles. 


Example 3 is a calculation in the deformed Ginibre of type considered in 27[, a model inspired by the 
doubly-degenerate Wishart ensemble We compute an integral representation and compare to known 

results in the vanishing external source limit. 

Example 4 serves as a proof-of-concept in applying the method to multiplication of matrices which 
attracted a lot of attention recently [s^ ■ In this case we calculate an integral representation for the averaged 
characteristic polynomial. 

Last example is a toy-model used to study the crossover between real /3 = 1 and complex /? = 2 Ginibre 
ensemble inspiried by the elliptic ensemble modelling the GUE-Ginibre transition. We arrive at the large N 
formula of the real-axis bump developed as we vary the crossover parameter. 


4.I. Ratio of determinants for complex /3 = 2 hermitian matrices 

In this example we show how to obtain formulas for averaged ratio of determinants by diffusion meethod. 
First, we introduce diffusion of the form 


drP{H,T) 




( 12 ) 
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where by convention we have Hij = + iyij) for i ^ j and Ha = Xu. We consider the ratio of 

characteristic polynomials: 


D{z,w,H) = 


det{z — H) 


det(ic — H) 

The extension to D consists of two Grassmann parameters p, q forming a super-determinant: 


D (z, w,H;q,p) = sdet 


w — H 


det{z — H — p{w — H) ^q) 


P z — HJ det{w — H) 

To solve condition (|4|) we open the superdeterminant with anti- ( 77 ) and commuting part: 


(13) 


(14) 


D 


d^drjexp 


-(e v) 


w — H 


p z — HJ \rj 
= + Ci^i) - ^ViVi - - fjip^i - ^,qr]i) + 


= / d^d-qe^^, 


-L ^ ^ 

fo 'y ' Xijifjiqj — qifjj + + ^i^j) H ^ y ) VijiVidj + ViVj + ~ 


i<j 


V2 


i<j 


where the exact proportionality factor is not needed in the calculation so we omit it. We calculate the LHS 
of (01) as 


AffD = ^ / d^dq I “ ViVj){'n^VJ + C*0) 


.Tg 


and we infer the laplacian Aa as 


(15) 


i<j 


Aa = ^ {dyjw - - 25p5q). 


(16) 


The final diffusion equation is therefore equal 


5 - 7 -D-j- — 2 A^ ^zz D-j-( 2 ;, w ^ q ') . 


(17) 


We comment on two peculiarities of (HZl). In the z direction it has a negative diffusivity constant and the 
diffusion happens also in the p,q Grassmann ’’directions”. Although the negative diffusion constant is not 
very physical, in the RMT context it is known to be responsible for universal oscillatory type behaviour Q . 
To deal with it on a technical level we simply Wick-rotate the z variable. As to the Grassmann diffusion we 
make use of the ’’flatness” property = 0,p^ = 0 and expand D(z,ty, i7;p, g) = + pG^^'^ -|- gD^^^ -|- 

gpD^^\ Using this decomposition we rewrite (ITT)) as an equivalent system of 4 equations for each 


O^IDt - 2A^ \^WW ^ZzJ^^T ’’ 


(18) 

(19) 

( 20 ) 
( 21 ) 


To find the solution of (HZl we first observe that only equations (flSl) and (I^Tl) contain relevant components 
i = 1,4. This is evident by observing that ultimately we are interested in the limit limp_q_>o D = We 
use a kernel of the heat equation 


N 


2ttt 


N 


N 


Kriz,w,y,v) = -—exp - — {v-wy - —(y-iz) 


2t 


2t 


( 22 ) 
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The solution to (EU is 


D 


-( 4 ) 


{z,w) = J dydvKr{z,w;y,v)T>''Q\-iy,v-Ho)=: (^Kr {z,w), 


(23) 


with Hq denoting the initial matrix and introducing the convolution operator o. With this notation, the 
solution to the inhomogeneous heat equation (fTSl) is 


= (^Kr O {z,w) = {z,w) - ^lDr^^\z,w). 

To write explicitly the solution we expand the initial condition 


D = 


det( 2 ; — Hq) 
det(r(; — Hq) 


1 + gpTr 


{z - Hq){w - Hq) J ’ 


n^^\y,v,HQ) = 


det(y - Hq) 


t^( 4)^ tr ^ det(y-iJo)^ 
Do '{y,v;Ho) = -r^Tr- 


det(n-i?o)’ ’ ’ det(z;-iJo) {y - Hq)(v - Hq)' 

We set diagonal matrix Hq = diag(/ii, h^)-, h = (hi,..., h]\j) and introduce 


44=n(-“ - A.) 


df^iw) = 


N 

2'kt 


N 




{q-hi)' 


The averaged ratio of characteristic polynomials (ETl) is equal to 


Dr{z,w) = TTj^{z)0j;iw) - ^J2^h-eS^)^h+e,M 


N 


2=1 


(24) 


(25) 

(26) 


(27) 


with Ci = (0...0,1,0, ..0) denoting a unit iV-dimensional vector in the i-th direction. Using (flUl) . we write 
the average explicitly as 


^{z,w)=C f 

^ ^ J det{w-H) 


(28) 


where tt and 9 are 


These types of averages are present as building blocks of biorthogonal structures [31 , 
the multiple orthogonal polynomials of type I and II respectively. 

To recover known formulas for the GUE case, we set hi = 0. Then 9k(w) = ^k-ifk-i(w) is proportional 

to the Cauchy transform fk^z) = J 
rewrite as: 


7rfe(s) and 7fe = n (f') V sf?- Along with = y, we 


Dr{z,w) = ^N-inN{z)fN-i{w) - t^n fN{w)TrN-i{z), (29) 

which is the ratio formula calculated for GUE ensemble in (l^ . 

4-2. Duality formula for complex (3 = 2 non-hermitian matrices 

Let Xij = Xij + iyij be a complex N x N non-hermitian matrix. We introduce diffusive dynamics of the 
form 

drPiX,r) = j^TrdxxfP{X,r) = ^X.t). (30) 

'i’J 
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We aim at calculating an averaged product of characteristic polynomials 


(Z, X) =det 


l[{zi-X)izi-X^) 


for which a deformation is a 2kN block matrix of the form 

0 Ipf — Ifc 0 X 


D('=)(Z,X;X) = det 


A0 1 


N 


-A^ 0 In 
Zl 0 Iat - Ifc 0X1 


( 31 ) 


(32) 


where Z = diag( 2 :i, ...Zk) and A is a complex k x k matrix representing the A-parameter space. We baptise 
the k-extended averaged characteristic polynomial (k-EACP) in agreement with [1^ where authors 
considered particular case of fc = 1. In the limit A —0 we recover the 
To proceed further, we open the using Grassmann variables 

D('=)(Z,A,X) = J dr]d^ e'^^, 

tg = Y. +E - E ^ - E - E +E A 


2=1 


2=1 2=1 
and 

'^xy — '~^xy H” '^^xy' 


2=1 






with kN dimensional vectors and {% = j = 1...N) and ” • ” denotes a sum over N dimensional 
indices. We set also A^v = a^v + We find that 


^ y drjd( I E 


v*j=i 


which determines the parameter laplacian as 

1 


Aa = ^Ti-a^At = ^ E {^Ij + 4) ■ 


(33) 


(34) 


We thus arrive at the hnal equation for the average = Dt-^^^(Z, A): 


9.d/''\z,A) = lTr9^^tD/''^(Z,A), 


(35) 


where we observe a reduction in dimensionality to k x k diffusion. The solution to this equation is written 
explicitly as 


D^^'''\z,Xo;A) = (^4^ J dBe-^'^dB-A){B^-A^)j^{k)^2:,Xo;B), 


(36) 


where A and Xq are the initial values of parameter- and the randomized matrix respectively. We turn to 
the product of characteristic polynomials by setting A —>■ 0: 


D. 


-{k), 


(37) 


’(Z,Xo)= (^4^ y’dSe-^™'DW(Z,Xo;B). 

To arrive at the duality formula, we use the explicit jPDF formula (HH) on the average D 

'^^''\z,Xo) =C' J dXe-^'^^^-^°'>^^^-^°'>D^'^\z,X). (38) 


-(fc) 
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We can therefore write the duality from (1571) and ((551) : 

C J + Xo;A = 0)= J dBe-^'^^^^'D^^Hz,Xo]B), (39) 

with the definition repeated for clarity: 


= det 


lAf — Ifc ' 
A Z) Itv 




zt 


-A^ (g) Iat 
g) Iat — Ifc (g 


xt 


(40) 


This new result is an extension of a similar formula for Xq = 0 obtained in M- Such duality type formulas 
were studied in hermitian ensembles by 241, for general /? in and in the context of string theory by 
among others. 


4-3. Deformed complex /? = 2 non-hermitian matrix 

In another example we deal with a variance deformed version of non-hermitian matrix model. Consider 
a standard diffusion equation 


a.P(X, r) = ^Y. (41) 

We are however interested in a slightly modified determinant of the form 

{z, X) = det {z - r-^Xfl-i) det {z - D-^X^T-^) . (42) 

This peculiar definition introduces a deformation of the usual non-hermitian model inspired by Wishart 
ensembles with time (T) and space (fl) correlations [2^. It arises by changing X —rXfl in both (l4lT) and 
dH 


P{X, t) ~ P(rxn, r) = exp 


|^_^Trr2(X - Xo)02(xt - 4)^ , 


(43) 


where Xq = T ^XqH With such reparametrization the characteristic polynomials (HU) averaged over 
standard non-hermitian diffusion (1421) becomes a standard characteristic polynomials averaged over doubly- 
correlated diffusion of type (1431) . To retain probabilistic interpretation, correlation matrices T and fl are 
positive definite and can be set diagonal. The extended determinant reads 


D(^’^)(z,X;r(;) = det 


z-T-^xn-^ 

n-^w 


-T-^w 


(44) 


with lim^_j.o We open it as 


Tq = r] ■ rjz + ^ ■ ^z — r] ■ {r ^Xfl ^) ■ rj — ^ ■ D ^X^T ^ ^ — rcTy • T ^ • ^-|-• fl ^-77, 

with N dimensional Grassmann vectors fj and rjj (j = I...X). Let the laplacian act on the deformed 
determinant 
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by condition (U) we infer that 




( 46 ) 


The final equation for the averaged (D)^ = reads 

drDr(z,A) = ^d^ii,Dr(z,A). ( 47 ) 

The equation is two dimensional so the solution reads, after w —>• 0 

_ IV 

D^(z,Xo) = — 

TTT 

D(^’^^(z,X;uj) = det 

which is the sought integral representation valid for general Xq and correlations F, 17. This studied average 
is given, using dm, explicitly as 



( 48 ) 


—r ‘^w \ 

( 49 ) 

Tl~'^w 

z-o-ixtr-d ■ 


Dr{z,Xo)=C' J dX exp (^-^Ty{X - Xo)HX - Xo)^ D^^’^'>{z,X), 
In the special Xq —?► 0 limit, the solution (HSl) reproduces the result of [27|: 

= — [ dppe-^P" 


(50) 


(51) 


Two complex Ginibre multiplied 

In this example we show that the method is suitable for product of matrices which recently had drawn 
much attention recently [^. We introduce two matrices undergoing an entrywise diffusion 


drP = ^Tra^ ^^P+ ^Tra^ ^tP, 

TV ^1^1 TV ’ 

and consider a determinant 

D{z,Xi,X2)=dei{z-XiX2)det{z-{XiX2)^) . 
To proceed with the method, the extension of D is an object linearized in W’s: 


D(z, Xi, X 2 ; u, V, w) = det 


(52) 


(53) 


( z 

—w 

0 

Xi\ 

V 

z 


0 

0 

xl 

u 

w 

U 2 

0 

— V 

u J 


(54) 


where one can check that lim„_>i_^^u,_>o D = D. We open the 4iV deformed determinants with the help of 
and t] 2 - 

D = J d^dpe^° , 

Tg = Ii- + 6 • ^2Z + fji ■ Piu + fj2 ■ r]2U - 6 • C2W + 6 • + m ■ r]2W - fj2 ■ r]iv+ 

+ ■ Xi ■ 1J2 + ^2 ■ xl ■ Pi + fji ■ xl ■ ^2 + 112 ■ X2 ■ ^ 1 - 
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The laplacian acting on D is 


^Xi,X2 = J^j (^1 ■ ^2»72 -rii+m- m^i ■ 6) e^'^, 
which dictates the parameter laplace operator to be of the form 

^ {dw,W “t“ . 

The final equation for averaged determinant reads 

dr'Dr{z;u,v,w) = ^ {dui^y, + dy^y) 'Dr{z\u,v,w), 


We write down the solution for v,w ^ 0 and rt —> 1: 

2 » 


Dt {z, (-’fi)o, (-’^ 2 ) 0 ) — 
with the definition 


N 

■nT 


'D{z, Xi,X 2 ; u, V, w) = det 


-f(l^ 


D(2 

,iXi) 


—w 

0 

XA 

V 

z 


0 

0 

xl 

U 

W 

VX2 

0 

— V 

U j 


(55) 


(56) 


(57) 


(58) 


(59) 


As before, we investigate the vanishing source limit (^i)o 0 where 

D(2:, 0, 0; w, n, ni) = (1 + vw)^ (|zp + vw)^ . 

The angles oi w,v in (l5^ can be integrated out with the help of hypergeometric function 


WAz) = dpdq gpe-"(«^+^>^)(|zp + 9V)"^2 Fi ( 


1 -Af _iv 1 4 |z|Vg^ 


2 ’ 2 ’ ’ (|z|2+p2g2)2 J ’ 

and we simplified further by introducing = ta, A = and integrating over a’s: 

—, , (2NV , f2Nt\ ,2 2^v fl-N N 4|zPt2 \ 

I —)(l = l +‘) 2 F.( — (|,((t.) 2 ). 

where the average, according to (ED, is explicitly given as 

'D;{z) = {C'f J dXidX 2 exp ^-^Tr(x|Xi + xlX 2 )^ det{z - X 1 X 2 ) det(z - (ATi^a)^). 
To the best of our knowledge, this result was not considered previously. 


(60) 


(61) 


(62) 


(63) 


4-5. Real/complex non-hermitian crossover model 

The diffusive method is applicable also to /3 = 1 real non-hermitian ensembles, we use it to show a simple 
crossover model between real and complex non-hermitian ensembles. An entrywise diffusion is defined by 


drP{X,T) 


I ^ f 2 \ 

dylm) 


P{X,t). 


(64) 


10 











By varying a between 0 1 the transition between /3 = 1 and /3 = 2 happens. We are interested in 

characteristic polynomial in the standard form 


D(z) = det(z — X) det(z — X^), (65) 

with N X N matrix Xij = xij + iyij. The extension is given by 

which, using the same techniques as previously, we arrive at an equation for the average characteristic 
polynomial 


_ 1 -|- _ 

-y (^Z ^ w'j (z, "Ui) . 


2N 


(67) 


The solution, after taking the w —>■ 0 limit, reads 


2N 


Dt(z,Xo) =- / drre D(z,r). 


( 68 ) 


which is valid for any external source Xq. For vanishing external source Xq —> 0, the formula (1681) agrees 
with results for real and complex Ginibre ensemble [3,81 Il5| . 

We are interested in the crossover region as it develops a bump located near the real axis z = 
as we also scale the crossover parameter a = 






2 

V2J’ 


(69) 


centered around 77 = 0. This deformation is interpreted as the precursor of a division between real and 
complex eigenvalues present in pure (3 = 1 Ginibre ensemble. 


5. Conclusions 

The presented method is shown to be universally applicable to general gaussian type random matrix 
models and serves as a calculational tool for obtaining ratios and products of characteristic polynomials. 

Examples considered present many new results of integral representations. In particular, a novel duality 
formula for products of characteristic polynomials in complex Ginibre matrices, a previously not considered 
characteristic polynomial for the product of Ginibre matrices and a real/complex Ginibre crossover model 
was considered. 

Main advantage of the method is a huge reduction in the degrees of freedom. It is applicable also 
to instances of statistical models with non-trivial covariance structures. Alongside, external source-type 
perturbations are present in this approach which broadens the general applicability of the method. 
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